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ABSTRACT 

We compute the virial mass of the central black hole (M) and the luminosity-to-mass 
(L/M) ratio of « 300 low-z quasars and luminous Seyfert 1 nuclei. We analyze: (1) 
whether radio-quiet and radio-loud objects show systematic differences in terms of M 
and L/M; (2) the influence of M and L/M on the shape of the H/3 broad component 
line profile; (3) the significance of the so-called "blue outliers" i.e., sources showing 
a significant blueshift of the [Om] AA4959,5007 lines with respect to the narrow com- 
ponent of H/3 which is used as an estimator of the quasar reference frame. We show 
that M and L/M distributions for RQ and RL sources are likely different for samples 
matched in luminosity and redshift, in the sense that radio-quiet sources have system- 
atically smaller masses and larger L/M. However, the L/M ratio distributions become 
indistinguishable if 8.5< logM< 9.5. Line profile comparisons for median spectra com- 
puted over narrow ranges of M and L/M indicate that a Lorentz function provides 
a better fit for higher L/M sources and a double Gaussian for lower L/M values. A 
second (redshifted) Gaussian component at low L/M appears as a red asymmetry fre- 
quently observed in radio-loud and radio-quiet sources with broader (FWHM>4000 
km s _1 ) H/3 broad component profiles. This component becomes stronger in larger 
mass and lower L/M sources. No specific influence of radio loudness on the H/3 broad 
component profile is detected, although equivalent widths of H/3 broad component and 
especially of [Om] AA4959,5007 are larger for radio-loud sources. We identify five more 
"blue outlier" sources. Since these sources are, on average, one magnitude brighter 
than other AGNs with similar mass, they are accreting at an Eddington ratio that is 
2-3 times higher. We hint at evolutionary effects that explain some of these results 
and reinforce the "Eigenvector 1" correlations. 
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1 INTRODUCTION 

Studies of emission lines play an important role in our under- 
standing of Active Galactic Nuclei (AGN) and the physics 
the accretion processes, although a connection between the 
main theoretical parameters (i.e., black hole mass, accretion 
rates, black hole spin, and a viewing angle) and observed 
spectral parameters is, at best, just sketchy and qualita- 
tive. In recent years, a major innovation has been ascribed 
to the so-called "Eigenvector 1" (El) correlations (Boroson 
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& Green 1992). The original El was most closely related 
to the anti-correlation between Fe n op t strength and peak 
[Om]A5007 intensity, and FWHM of the H/3 broad compo- 
nent (H/fec). 

El-related correlations have emerged as very robust, 
since they continue to appear in different low- 2: AGN sam- 
ples, even if a set of spectral parameters partly different 
from the ones used for the original Principal Component 
Analysis (PCA) by Boroson & Green (1992) is considered. 
An outstanding relationship in the context of El involves 
the FWHM(H/3bc) and the ratio between the equivalent 
width of the Fellopt complex centered at A4570 and H/3bc, 
Rpoii= W(FeiiA4570)/W(H/3 B c). Low redshift AGN occupy 
an elbow sequence in the plane defined by these two pa- 
rameters (see Sulentic et al. 2000a; Marziani et al. 2001). 
The overall sequence is defined by RQ sources while the 
RL AGN occupy only a restricted part of it. The domain 
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where FWHM(H/3 B c)< 4000 km s _1 is predominantly oc- 
cupied by RQ sources (Sulentic et al. 2000a; see also Su- 
lentic et al. 2003). There is an apparent dichotomy in the 
H/3bc line profile shapes: sources with FWHM(H/3bc) ;$ 
4000 km s _1 show typically a Lorentzian H/3bc, while for 
FWHM(H/3bc) > 4000 km s" 1 a double Gaussian fit is 
more appropriate (Sulentic et al. 2002); (3) the so-called 
"blue outliers" (sources showing a significant blueshift of 
the [Om] AA4959,5007 lines with respect to inferred quasar 
frame) occur among sources with strongest FellA4570 nar- 
rowest H/3bc (Zamanov et al. 2002; Marziani et al. 2003a). 

The robustness of the El correlations stems probably 
from a series of concurring factors: (1) an Eddington ratio 
that may be continuously changing along the El sequence 
(e.g., Boroson, Persson & Oke 1985; Boroson & Green 1992; 
Murray et al. 1995; Nicastro 2000; Marziani et al. 2001; 
Boroson 2002); (2) outflow/wind properties which are ex- 
pected to depend on Eddington ratio (e.g., Dultzin-Hacyan 
et al. 2000; Marziani et al. 2001); (3) orientation effects that 
will blur any physical correlation and displace sources in 
the direction of increasing Eddington ratio (Marziani et al. 
2001; see also, Nagao et al. 2000, Zhou et al. 2003); (4) 
Broad Line Region (BLR) structural effects (Sulentic et al. 
2000b, Sulentic et al. 2002; see also the Proceeding of the 
Nebraska meeting for the El relevance for the Broad Line 
Region structure; Gaskell et al. 1999); (5) last, but not least, 
evolutionary effects that will lead to an increase of black 
hole mass as well as of the [Olll]A5007 rest-frame equivalent 
width (Zamanov et al. 2002; Boroson 2002). We will briefly 
discuss them in § 18.41 and § 9. 

Even if there is not yet a consensus (see e.g., Marziani et 
al. 2001; Boroson 2002), the physical drivers of El are likely 
to involve the luminosity-to-black hole mass ratio (L/M oc 
Eddington ratio), an orientation angle, and the black hole 
mass (M) (Zamanov & Marziani 2002). Additional proper- 
ties such as black hole spin (e.g., Wilson & Colbert 1995) or 
the properties of the host galaxy environment (e.g., McLure 
et al. 1999) are also likely to play a role in any RQ-RL di- 
chotomy. 

This paper presents a basic investigation of the influ- 
ence of L/M and M on RL and RQ sample properties, on 
the H/3bc line profile, and on the occurrence of the "blue 
outliers", which are all major El-dependent phenomena. In 
§ 2 we present our quasar sample followed by (§ 3) derivation 
of M and L/M, (§ 4) presentation and discussion of mass - 
luminosity diagrams, (§ 5) a comparison of radio quiet and 
loud sources, (§ 6) an analysis of H/3bc profiles, (§ 7) as well 
as significance of the "blue outlier" sources. Our results are 
discussed in § 8. § 9 provides some concluding remarks. 



2 SAMPLE DEFINITION & MEASUREMENTS 

We use spectra of 278 AGNs. Our AGN sample includes 
215 sources and has been presented elsewhere (Marziani et 
al. 2003b, hereafter M03). We added 63 soft X-ray selected 
sources with suitable optical spectroscopic data (Grupe et 
al. 1999; hereafter G99) because this turns out to be the 
most fertile source of "blue outliers". Inclusion of Grupe et 
al. data makes our sample less biased, as described in H5.4I 
The optical characteristics of these X-ray selected sources 
include: (1) relatively "narrow" broad emission lines (e.g. 
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Figure 1. Distribution of absolute B magnitude (Mb) versus 
black hole mass M (in solar units) for the 278 AGN sample. Mb 
values are corrected for the galactic extinction. Triangles and cir- 
cles indicate radio-quiet and radio-loud sources respectively. Open 
symbols indicate sources with FWHM(H/3bc)^4000 km a -1 ; 
filled symbols sources with FWHM(H/3 B c)<4000 km s _1 (see 
85.41 . The bars in the lower right corner indicate typical errors. 

Narrow Line Seyfert 1 [NLSyl] sources), (2) strong optical 
Fe ilopt blends and (3) weak forbidden emission lines. The 
merged datasets provide reasonable quality spectra (S/N~ 
20 - 70) for the region near H/3 in 278 AGN (z < 0.8). 
The M03 dataset is not complete (see M03 for discussion of 
completeness and sample biases). 

We use line parameters reported in M03 and de- 
rive the corresponding measures (using the same pro- 
cedures described in M03) for the X-ray selected sam- 
ple from the digital spectra made publicly available by 
G99. We consider here measurements of the following pa- 
rameters: (a) FWHM(H/3bc); (b) equivalent width ratio 
Rpoii=W(FeiiA4570)/W(H/3 B c), (c) FWHM measurement 
of the lines in the FellA4570 blend (see M03 for details), 
(d) radial velocity difference between [Om]A5007 and H/3 
at the line peak (without narrow component subtraction). 
Accuracy of the G99 wavelength calibration was evaluated 
by examining the distribution of radial velocity difference 
between [OIII]A4959 and [OIIIJA5007. Whenever possible, 
[OIIIJA4363 was also used and all inter-narrow line differ- 
ences were found to be comparable to our measures (cf. Za- 
manov et al. 2002). 

Apparent (my), and absolute B magnitudes (Mb) were 
taken from 10 th edition of the AGN catalog (Veron-Cetty 
& Veron 2001). A magnitude correction +0.88 mag was ap- 
plied to convert Mb from Ho=50 km s _1 Mpc~* to Ho=75 
km s _1 Mpc -1 (qo = 0) which we use throughout this pa- 
per. Galactic extinction values Ab were taken from NED 
following Schlegel et al. (1998). Sources were considered RL 
if the specific flux at 6 cm and in the B band was larger than 
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10, with data taken from the Veron-Cetty & Veron (2001) 
catalog. 



3 COMPUTATION OF M AND L/M 

3.1 Black Hole Masses 

One can estimate the mass of the putative supermassive 
black hole using FWHM(H/3bc) and a reverberation "ra- 
dius" (Kaspi et al. 2000) along with the assumption of viri- 
alized motions. It is now common to estimate the mass 
by assuming the BLR distance from the central continuum 
source tblr oc (Lsioo)" and a =0.7, as derived from the 
reverberation data (Kaspi et al. 2000). The virial mass is 
M=r B LRV 2 /G, where v= %/3/2 FWHM(H/3 BC ) (c.f. Woo & 
Urry 2002a), and G is the gravitational constant. We there- 
fore can write the black hole virial mass as follows: 

M = 4.817 x ( ™M(H/M W A L 5100 V 7 

\ lkms- 1 J \W 44 ergs- 1 J ' w 

where M is in solar units, L5100 is the specific luminosity at 
5100 A (in units of ergs s _1 cm -2 A -1 ), and 

AL5100 = 3.137 x io 35 -°- 4 ( m b-Ab) ergg -i (2) 

In order to convert the absolute B magnitude Mb to AL5100 
we assumed that the specific flux is f„ oc u~ 3 . The derived 
masses are thought to be reasonable estimates within a fac- 
tor of about 2-3 (see also Woo & Urry 2002a, Vestergaard 
2002). 

3.2 Influence of Orientation 

Source viewing angle is especially important in AGN where 
the concepts of disk accretion/line emission and jet ejection 
are widely favored. In the above mass derivation formula we 
did not include any orientation related broadening factor. If 
H/3bc is emitted in a highly flattened structure then M can 
be underestimated by a factor ~ sini" 2 in pole-on/face-on 
sources (i.e., with accretion disk axis oriented along our line 
of sight). Unfortunately, we do not yet have a proper way to 
estimate the viewing angle i for each source. This is a major 
drawback because there is observational evidence suggest- 
ing that orientation can affect line width by a factor ~ 2 
(Brotherton, 1996; Marziani et al. 2001; Jarvis & McLure, 
2002). Even though we neglected the effects of orientation in 
the general sample population, we considered it in compu- 
tations of M and L/M the relatively rare and extreme class 
that is likely observed close to "pole-on" : the "blue outliers" 

(JH- 

3.3 The Role of the BLR Size-Luminosity 
Relationship 

The exponent in the relationship between source luminosity 
L and tblr is assumed to be a=0.7. Any deviation from 
this value has quantitative effects on our results. Our sample 
includes many moderately luminous quasars at 0.4 < z < 
0.8, so it is important to stress that the relationship derived 
by Kaspi et al. (2000) is based exclusively on quasars of 
z < 0.4 and that the high (and low) luminosity ranges of 
the correlation are poorly sampled. 



If we consider sources in the luminosity range 43.4 < 
logL/L0 < 45 (i.e. where we have uniform luminosity sam- 
pling), the slope of the best fit is q=0.8, and could be eas- 
ily as high as a = 1 without increasing significantly the fit 
standard deviation. Marziani et al. (2001) considered this 
case which seems appropriate for the PG quasar luminos- 
ity range. One must remain open to the possibility that 
0-5 i$ ol < 1, and that a might even be a function of L. 
Changing a implies an L-dependent change in mass esti- 
mates; therefore the slope of the luminosity-to-mass rela- 
tionship (Fig.0 is affected as well as the location of points in 
the L/M vs. M diagram (Fig.[7|& Fig.E^. Despite these pos- 
sibilities, systematic trends discussed in this paper should 
not be affected. 

3.4 Bolometric Luminosity 

We calculated the bolometric luminosity L from L ~ 
10A L A (5100A) (details can be found in Wandel, Peterson, 
& Malkan, 1999, Elvis et al. 1994, Collin et al. 2002). Values 
computed in this way were used to derive the L/M ratio. The 
L/M ratio everywhere is expressed in solar units with the 
solar value (L/M)© = 1.92 ergs s -1 g" 1 , and the Eddington 
limit corresponding to logL/M = 4.53. 

We neglected any possible differences in the spectral 
energy distribution between source subclasses. In order to 
test this benign neglect we considered 44 sources that are 
common with Woo and Urry (2002a) where SEDs were con- 
structed from archival data. Comparison of our bolometric 
luminosity values showed no significant systematic differ- 
ence: AlogL « —0.06 for all 44 sources, with a standard 
deviation a « 0.24. RQ sources show A log L ~ —0.06 with 
a « 0.22, while RL AlogL « -0.07 with a slightly larger 
a ~ 0.26. If we consider separately RQ and RL sources with 
the restrictio n FWHM(Hfec)^ 4000 k m s" 1 (see fOl . we 
find AlogL w -0.05 (RQ) and AlogL w -0.01 (RL). It is 
worth noting that the scatter is due to a minority of bad- 
behaving data points («20%). If they are removed, the stan- 
dard deviation becomes a ~ 0.1 in all cases, with systematic 
differences always <0.05. 



4 MASS - LUMINOSITY DIAGRAM 

FigE shows a plot of Mb vs. M for our combined sam- 
ple which covers an absolute B magnitude range —20 < 
Mb< —27 and an estimated black hole mass range 
M < 10 (M is everywhere given in Mq ). The data show 
a rather well-defined range in L/M with almost all sources 
lying between 0.02 < L/LEdd < 1.00. RQ sources show evi- 
dence for significant Malmquist bias while RL sources show 
the opposite trend probably related to a bias towards se- 
lecting higher luminosity core-dominated sources likely to 
be beamed. 



5 RADIO LOUD AND RADIO QUIET AGN 

Table^presents average and median values of relevant prop- 
erties for the RL and RQ subsamples. Kolmogorov-Smirnov 
tests indicate that RL and RQ samples differ significantly 
in redshift, my and absolute B magnitude distribution. The 
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Table 1. RQ/RL sample averages and medians for redshift, apparent V magnitude (my), absolute B magnitude (Mb), BH mass (M), and 
luminosity-to-mass ratio. Average values are given with sample standard deviations, median values with first and third quartile values. 
Pks i s the probability that the two distributions are randomly drawn from the same parent population using Kolmogorov-Smirnov tests. 
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Figure 2. Distribution of apparent V magnitude (left panels) 
and redshift z (right panels) for our 278 AGN sample. Shaded: 
RL; Unshaded: RQ. The left panel shows that, in our sample, RL 
sources are on average less bright than RQ sources while the right 
panel suggests a distribution favoring higher z for RL AGN. 



distribution of the computed parameters M and L/M are 
also significantly different for RQ and RL samples. 

RL sources are somewhat over-represented in the M03 
sample, but not in the M03+G99 sample. Any over- 
representation can be properly quantified if we compare gen- 
eral population expectation for the number ratio between 
the number of RL sources and the total number of sources 
(fn =N(RL)/[N(RL)+N(RQ)]) in the M B range covered by 
our sample. Weighting the luminosity-dependent f R value 
from La Franca et al. (1994) over the Mb distribution of our 
sample (see Fig. 2 of M03) i.e., within -27 < M B < -23, 
we obtain f R « 25% versus f R « 39% in M03 and f R « 26% 
in M03+G99. 

More cumbersome are the biases shown in Fig. [5] RL 
sources tend to be fainter and more distant than RQ ones, as 
well as intrinsically more luminous. Biases that are affecting 
RL sources are not clear. RL core-dominated sources are 
over-represented in M03. This will contribute to bias the 
RL sample to higher L and z values (see Sulentic et al. 2003 
for a more thorough discussion). 
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Figure 3. Distribution of M (left panels) and L/M ratio (right 
panels) for our sample (upper panels) and for the sample by 
Grupe et al. 1999 only. Shaded: RL; Unshaded: RQ. The left 
panel shows that RL sources shows larger M than RQ while the 
right panel suggests a largely overlapping L/M range for RQ and 
RL AGN, even if the largest L/M sources are almost all RQ. 



5.1 Bootstrap Simulations 

Are differences in redshift and luminosity distributions the 
cause of the inferred M and L/M differences between RQ and 
RL sources? We can use most of our sources in creating bet- 
ter matched subsamples out of our M03+G99 sample (boot- 
strap techniques, i.e. Efron & Tibshirani, 1993). The most 
important parameter appears to be redshift (which implies a 
strong bound on the apparent magnitude). We constructed 
pseudo-samples of RQ and RL sources with similar redshift 
and apparent magnitude distributions. More precisely: (1) 
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of the order of unceratainty in M and L/M determination 
and therefore too small to be of relevance). The similarity 
of the parameter space occupation for RQ and RL sources 
with FWHM(H/fec) > 4000 km s" 1 carries into a similarity 
in derived L/M and M. In other words, there is a large RQ 
population whose M and L/M values are similar to the ones 
of RL sources. 



Figure 4. Results of bootstrap simulations. Distributions of av- 
erage M (left) and L/M (right) for 1000 pseudo samples of 40 RL 
(shaded) and RQ objects with matching apparent magnitude and 
redshift distribution. Bootstrap simulations suggest that there are 
systematic differences in the sense that RL sources have large M 
and lower L/M ratios than RQ AGN. 

we randomly selected two pseudo-samples with the same 
z and my distributions, indistinguishable within a 2a con- 
fidence limit, (2) we computed average logM and logL/M 
values for each pseudo-sample, (3) we repeated the selection 
~1000 times, (4) we computed the distribution of average 
logM and log L/M values for the pseudo-samples (the boot- 
strapped samples), (5) we d erived the expectation values for 

< logM > and < logL/M > as the medians of the distri- 
butions for the bootstrapped samples, (6) we estimated the 
significance level of any difference from the average values 
and the dispersion of the distributions. While cannot elimi- 
nate biases from our sample, we can try to apply the same 
bias to both populations. 

5.2 Matching z and niv Distributions 

The bootstrap procedure makes both the M and L/M dif- 
ferences stronger (see the upper panels of Fig[3]and Fig|3J: 
the median values are < log M(RQ) > ^ 8.07 ± 0.07 vs. 

< logM(RL) >w 8.81 ± 0.06, and < logL/M(RQ) >^ 
3.88 ± 0.05 vs. < log L/M (RL) >^ 3.45 ± 0.06. The results 
can be seen in Fig|3] (the original sample) and Fig^] (after 
bootstrap simulations). Since matching mv and z is equiv- 
alent to a match in L, it follows that the two plots are not 
independent, and that systematically larger FWHM(H/3bc) 
drives a larger mass estimate for RL sources. 

5.3 Same z and my Distribution with a Narrow 
Mass Range 

To avoid sample biases, and to assess the reality of any L/M 
systematic difference, we restricted our attention to sources 
within a narrow mass range. At the same time, we retained 
the above conditions on the z and my distribution. The 
mass range was chosen so that we had the maximum possible 
number of sources: 8.5 < log M < 9.5 with 48 sources (Fig. 
0. No systematic difference in logL/M (< 0.2) is found. 
In a comparison of RQ and RL sources with the additional 
condition FWHM(H/3bc) > 4000 km s _ \ and matching the 
z and my distributions, we find again no systematic dif- 
ference in either M or L/M distributions (see Fig. |S| note 
that the difference can be statistically significant, but it is 



5.4 Population A and B: Understanding Biases 

Sulentic et al. (2000a) introduced the concept of Popula- 
tion A [FWHM(H/3 B c) < 4000 km s" 1 ] and Population B 
[FWHM(H/3bc)>4000 km s" 1 ] as fundamentally related to 
the BLR structural properties (Sulentic et al. 2000a,b). RQ 
sources dominate Pop. A (88 % of the total) and are about 
| of Pop. B sources. Most RL sources («75%) are Popula- 
tion B. Thus to a large extent a RL vs. RQ comparison is 
a Pop. A RQ vs. Pop. B RL comparison (as it was in the 
early study on CivA1549 by Marziani et al. 1996). 

The validity of the previous results on M and L/M dif- 
ference depends on (1) the fraction of Pop. A RQ sources at 
low z, and on (2) whether most RL sources at z < 1 are Pop. 
B. A proper assessment of the Pop. A/Pop. B ratio for RQ 
samples would imply a thorough analysis of the discovery 
biases affecting major surveys (see e.g., Oshlak et al. 2002). 
This goes far beyond the aim of the present paper. 

Actually, we know that we have a bias favoring Pop. 
A RQ sources in the G99 sample. If we apply the same 
bootstrap analysis to the M03 sample only, we obtain a 
lower fraction of Pop. A sources (« 45%), but we reach 
the same qualitative conclusions. Therefore a significant bias 
may arise only if we miss a large fraction of Pop. B RQ 
sources (>20%). This does not seem the case. An analysis 
of the SDSS shows that NLSyls should be » 15% of all low- 
z AGN (Williams et al. 2002). The ratio between NLSyls 
and the rest of Pop. A (i.e., sources satisfying the criterion 
2000 km s~ 1 <FHWM(H/3bc)< 4000 km s" 1 ) is « § (es- 
timated on the whole M03 sample, on a restriction to PG 
sources). This implies that Pop. A should be « 60% of all 
low z AGN. We note that: (a) in M03+G99, Pop. A sources 
account for ~52% of all sources; (b) the ratio of NLSyls 
to all RQ sources is « 15-20% if we consider the M03 and 
M03+G99 independently, in line with the SDSS findings. 
Therefore, if there is a sample bias favoring one of the two 
populations, it is not strong. 

Regarding the second issue, it has been possible to 
show that Fanaroff-Riley II quasars belong exclusively 
to Pop. B. RL sources with narrower lines [RL with 
FWHM(H/3 B c)<4000 km s _1 ] are predominantly core- 
dominated (CD) sources, and therefore likely to be a less fre- 
quent, beamed population of preferentially aligned sources 
in a FRII/CD unification scenario (Sulentic et al. 2003). 

Even if the amplitude of the L/M and M systematic 
difference will of course depend on sample definition, such 
difference is unlikely to be caused by selection effects, and 
should be considered real. This is in line with recent research 
pointing toward systematic differences in the mass function 
of RQ and RL AGN (as discussed in flOl . 
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Figure 5. Distri bution o f AlogM = log M(RQ) - logM(RL) 
(left panel) and AlogL/M = log L/M(RQ) - log L/M(RL) (right 
panel) for ~1000 pseudo samples of RL and R.Q sources, my and 
2 distributions were matched, and a restriction to the M was ap- 
plied: 8.5 $C logMsC 9.5 (48 sources). RL and RQ quasars have 
almost the same L/M ratio distributions. 
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Figure 6. Distribution of AlogM (left panels) and Alog L/M 
(right panels) for ~1000 pseudo samples of 40 objects with 
FWHM(H/3 BC )^4000km s _1 but separated on the basis of radio 
loudness. It points that RL and RQ Pop. B are not distinguishable 
on the basis of M or L/M. 



6 H/?bc PROFILE SHAPE 

6.1 The Influence of M and L/M 

We generated median H/3bc profiles in bins covering narrow 
ranges of derived M and L/M. The adopted binning is shown 
in Fig. [7| The median spectra (normalized to the local con- 
tinuum) are plotted in Fig. |H| Table 2 provides, for every 
bin, average measures of the H/3bc profile centroid at 1/4 
and 1/2 peak intensity, the best fit to the profile shape, and 
the intensity ratio (C2/C1) in case when two components 
have been used. 

The median H/3bc profile computed for each bin (see 
Fig. Inland 11011 suggests that: 

• a Lorentz function provides good fits to median pro- 
files in bins for which 3.9 < log L/M < 4.4, and for 
3.4 < log L/M < 3.9 and 7.5 < logM < 8.5. Profile models 
employing Gaussian functions yield poorer fits with signifi- 
cantly larger \ 2 - The slight redward asymmetry in the bin 
3.9 < log L/M < 4.4 and 8.5 < log L/M < 9.5 can be mod- 
eled by adding a weak, redshifted Gaussian component C2 
with C2/C1 ~ 0.1. Similar considerations apply to the bin 
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Figure 7. Binning of our sample in the L/M vs. M parameter 
space for the median spectra computations. Symbols are as in 
Figure 1. 



4.4 < log L/M < 4.9 and 7.5 < logM < 8.5 (i.e., largest 
L/M, possibly super-Eddington, and lowest M) profile: a 
slight blueward asymmetry visible in the median profile can 
be modeled as an additional, blueshifted component con- 
tributing ~ 0.1 of the total line emission. 

• If log L/M < 3.4, profiles are better fit with double 
Gaussian models than using a Lorentz core + broader Gaus- 
sian component The latter model is favored if the narrow 
component of H/3 is not subtracted (e.g. Oshlack et al. 2002, 
but lack of narrow component subtraction will also lead to 
an underestimate of M. 

• The "red shelf or redshifted Gaussian component (we 
refer to it as very broad line region (VBLR) component) ap- 
pears to be strongly influenced by M, i.e. larger M sources 
show a more prominent VBLR component, as well as by 
L/M (this is especially evident from the values for 8.5< 
logM sC 9.5 reported in Table 2). The M dependence is eas- 
ily seen by comparing the solid (lower mass) and overlayed 
dashed (higher mass) profiles in each frame of Figure |H] Ta- 
ble 3 shows that c(l/4) measures are significantly redshifted 
for the bins corresponding to larger M. 

Best fits to H/3bc and individual line components are 
shown in Fig. [5] Fig. 1101 illustrates a case in which the ex- 
change from a Lorentzian to a Gaussian significantly worsen 
residuals. The case shown is the one where fits provides more 
similar results. In general exchanging functional forms from 
the best fitting ones imply larger differences and a \ 2 wors- 
ening by a factor of ~ 2. 



6.2 Influence of Radio Loudness on H/Jbc 

Is the H/3bc profile shape influenced by radio loudness? To 
answer this question we generated the H/3bc profile (see 
Fig. El fo r the interval 8.5 < log M < 9.5 and any L/M 
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Table 2. H/3bc li ne centroids at 1/4 and 1/2 fractional intensity as a function of logL/M. For each interval of logL/M the table reports 
the H/3bc best fitting function (FF; either Lorentzian (L) or Gaussian (G)), and their peak position and FWHM. In case two component 
are needed for the best fit, their intensity ratio (C2/C1) is reported in the last column. 
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for RQ and RL separately (Nrq = 56, Nrl = 36). Fig. 
1111 shows that the composite spectra of RL and RQ ob- 
jects are different: the narrow lines ([Olll]AA4959,5007 and 
H,0nc) are stronger and W(H/3bc) is lower in RL objects. 
At the same time, the H/3bc profiles remain indistinguish- 
able within our S/N limits (see Fig. llH . To further check this 
result we additionally generated the composite spectra sepa- 
rating RL and RQ objects in three other M and L/M ranges: 
(a) 8.5 < logM< 9.5 and 3.9< logL/M < 4.4 (Nrq =14, 
Nrl =19); (b) 8.5 < logM< 9.5 and 3.4< logL/M < 3.9 
(Nrq =26, Nrl =13); (c) and 8.0 < log M < 9.0 and 3.9 
< logL/M<4.4 (Nrq =14, Nrl =19). Although the S/N 
ratio of the composites was lower since fewer objects were 
used, results were similar: indistinguishable H/3bc profiles 
and centroids, with higher W([Olll]AA4959,5007) and lower 
W(H/3bc) for RL objects. We conclude that the H/3bc line 
profile shape does not depend strongly on radio loudness. 



7 THE "BLUE OUTLIERS" 

7.1 Spectral Properties of [Om] AA4959,5007 "Blue 
Outliers" 

In a previous investigation (Zamanov et al. 2002), we iden- 
tified seven objects from the M03 sample that showed an 
[Om]AA4959,5007 blueshift relative to the peak of H/3. The 
shift amplitudes were ^ 250 km s _1 leaving little doubt 
that we are observing a significant velocity displacement. 
We have identified five additional AGN that show large am- 
plitude blueshifts among the objects of G99. Some properties 
of all of these sources are given in Table [3] All of the sources 
show strong Fe n op t emission. In the El context they show: 
(1) RfcII>0.5, (2) FWHM(H/3) < 4000 km s -1 (i.e., they are 
located in spectral bins A2 and A3 of Sulentic et al. 2002); 
(3) a large ClvA1549 broad component blueshift (Zamanov 
et al. 2002). Blue outliers show also weak [Om] AA4959,5007 
emission with W([Olll]A5007)< 20 A, and a median value 
of «6 A. 



Nine out of our 12 outliers are formally NLSyls. In gen- 
eral, the [Oin]AA4959,5007 lines of NLSyls have a relatively 
narrow profile with often, in addition, a second broader, 
blueshifted component (Veron-Cetty, Veron & Gongalves 
2001). The blue-shift, is 400 - 1100 km s _1 when the two 
components can be deblended (Zheng et al. 2002). These 
values correspond to the blue outliers range suggesting that 
the latter are an extreme subset with low W([Om] A5007) 
when only the shifted component is visible. 



7.2 FWHM of Fell pt and [Om] AA4959,5007 lines 

Computation of M requires that H/3bc be separated from 
any narrow component so that FWHM(H/3bc) can be prop- 
erly measured. The problem is that extreme sources like 
the blue outliers show no profile inflection (reflected in the 
fact that they are well fit by a Lorentzian function). The 
rationale for our methodology has been detailed in M03, 
and will not be rediscussed here. We just remind that we 
normally assume (and find) that the narrow component 
of H/3 shows almost the same shift and width properties 
as the [Om] AA4959,5007 lines. However, in blue outliers 
W([Oiii]A5007) is smaller than in any broad emission line 
AGN class and is also blueshifted relative to the peak of H/3. 
We conclude that we are not measuring a narrow line com- 
ponent at the Lorentzian tip of the H/3 profile. Any narrow 
component of H/3 is assumed to be blueshifted and lost in 
the blue wing of the broad line (even assuming it has the 
same strength as [Oiii]A5007 in blue outliers it would rep- 
resent less than 10% of the total line flux) . Subtraction of a 
narrow H/3 component would have a non-negligible effect on 
the FWHM(H/3bc) measurement which is needed for the M 
estimation. 

The FellA4570 feature is a complex of lines thought 
to have FWHM similar to H/3bc- We estimated 
FWHM(FellA4570) independently from FWHM(H/3bc) 
by Fellopt template fitting using the I Zw 1 spectrum 
(of course, only for highest S/N spectra; see M03 for 
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Figure 11. Comparison of the RQ and RL objects in the mass interval 8.5 < logM< 9.5 and any L/M ratio (N^q =56, Nrl =36). 
(leff) Comparison of the composite spectra of RQ and RL (after subtraction of the Fe n op t template). W(H/3bc) is 20% higher in RQ, 
W([Oiii] A5007) is 35% higher in RL. (right) cleaned and scaled H/?bc profile. The H/3bc profiles of RQ and RL are almost identical. 



details). In our sample FWHM(FenA4570) is similar to 
FWHM(H/3bc) with few exceptions only (see Bongardo et 
al. 2002). Estimated FWHM(FeiiA4570) values are given 
in Table [3] No narrow FellA4570 component has been 
detected and nor is one expected given the low inferred 
density for the narrow line region. In principle, we can use 
FWHM(FellA4570) instead of FWHM(Hftc) in the mass 
derivations in order to avoid any error associated with a 
narrow line component. Use of FWHM(FellA4570) instead 
of FWHM(H/?bc) will yield M values within ~10% of ones 
derived from FWHM(H/? B c). 

The correlation of nuclear BH mass with stellar bulge 
velocity dispersion <r» is now well established in nearby 
galaxies. Recent results (Nelson 2000; Boroson 2002) show 
that a BH mass - FWHMQOm] A5007) correlation is also 
present but with larger scatter. The FWHM([Om] A5007) 
measures given in Table [3] provide another way to esti- 
mate BH masses. The derived values are however consid- 
erably higher than values calculated using FWHM(H/9bc) 
and source luminosity. This is not surprising because 
NLSyl-type AGN apparently do not follow Nelson's rela- 
tion (Mathur, Kuraszkiewicz & Czerny 2001). If NLSyls do 
not follow the relation then blue outliers will almost cer- 
tainly show the same lack of agreement. This does not con- 
tradict our assumption that blueshifted [Om] AA4959,5007 
arises in outflowing gas (Zamanov et al. 2002), possibly as- 
sociated with a disk wind. The [Om] AA4959,5007 region in 
blue outliers may be very compact and its velocity field is 
not likely to be dynamically related to the host galaxy stel- 
lar bulge. This points to a limiting W([Om]A5007) (w 20A) 
below which FWHM([Om] AA4959,5007) emission ceases to 
be a useful mass estimator. 



7.3 M B and L/M Ratio of [OIII] Outliers 

All blue outliers lie in the region of El thought to be pop- 
ulated by highly accreting (high L/M) sources (Marziani et 
al. 2001). Most of them are located close to the Edding- 



ton limit (log L/M » 4.53, see Fig[HJ. This is also true 
if: (1) FWHM(FenA4570) is used instead of FWHM(H/3 B c) 
to derive M, and (2) an orientation correction is applied. 
The blue outliers may be oriented close to pole-on. Since M 
oc FWHM 2 , an inclination correction will move the "blue 
outliers" toward higher mass and lower L/M (see the ar- 
rows on Fig |12ll . Even with an inclination correction by a 
factor ~4, the L/M values for blue outliers remain among 
the highest observed. We note that: (1) other sources should 
also be corrected for orientation especially if low ionization 
lines are emitted in a flattened configuration in all Popula- 
tion A sources, as suggested (Marziani et al. 1996) ; (2) some 
care is needed since we neglected a photometric correction 
which is basically unknown. We do not know whether any 
sort of continuum beaming operates for RQ AGN as i — > 0° 
(very similar W(H/3bc) for CD and LD RL sources suggests 
no strong effect; Sulentic et al. 2003). Observational errors 
and uncertainty in the orientation correction do not allow 
us to determine whether some blue outliers are really super- 
Eddington sources as they appear in Fig |12l 

All blue outliers show estimated masses less than 
1O 8,5 M0. We tried to determine if the blue outliers large 
L/M was a consequence of a low-mass bias by calculating 
the Mb difference between blue outliers and other sources 
in two mass ranges: (1) 7.15 ^ log M ^ 8.45 (where the 
"blue outliers" are before the orientation correction) and 
(2) 7.55 ^ log M ^ 8.85 (where they move after the orien- 
tation correction). Table [I] shows that the blue outliers are 
almost certainly ~ 1 mag brighter than other sources with 
similar mass implying 2-3 times higher L/M ratios. 



8 DISCUSSION 

8.1 Mass-Luminosity Diagram 

The ranges of M and L/M considered in this study are 
similar to those in other recent work (Woo & Urry 2002a; 
Collin et al. 2002). We find that the Eddington limit de- 
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Table 3. The spectral properties of the "blue outliers": radial velocity difference between top of H/3 and [Oin]A5007, equivalent width 
and FWHM of [Oiii]A5007, FWHM of FeoA4570, the equivalent widths ratio of FellA4570 and H/3 BC (R Fe Il)- 
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Table 4. Mb of the "blue outliers" (as they are given in Veron- 
Cetty &: Veron catalog without correction for the cosmology), 
compared with other objects in the sample with similar mass. Pks 
is the probability that the "blue outliers" and the other objects 
in the indicated mass ranges are randomly drawn from the same 
parent population according to Kolmogorov-Smirnov test. 
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fines an approximate upper boundary to the luminosity dis- 
tribution (as in Woo & Urry 2002a), indicating that there 
are no known low-z AGN accreting significantly above the 
Eddington limit (this depends on the adopted Ho value). 
Figure also suggests that there may be fewer high (than 
low) mass sources with L/M close to the Eddington limit. 
This might be a selection effect or an indication that galax- 
ies with a high mass BH may be unable to supply fuel at 
high L/M. In other words, the upper envelope in Figure 
might be due to the Eddington limit as well as the ability of 
the surrounding matter to feed the accretion flow (see, e.g., 
Nicastro et al. 2003). The lower envelope of the luminosity 
distribution may be due to a selection effect (e.g. Woo & 
Urry 2002a), or it may indicate that only sources radiat- 
ing at 0.01 < (L/LEdd) < 1-00 exhibit a stable broad lines 
region. The latter possibility might be connected with the 
presence/absence of an accretion disk wind (see also Kol- 
latschny & Bischoff 2002). 

8.2 On the Difference Between Radio Loud and 
Radio Quiet AGN 

Somewhat confusing claims have been recently made on the 
difference between RL and RQ sources in terms of M and 
L/M. Laor (2000) found that the radio loudness is strongly 
related to M. Woo & Urry (2002b) concluded that radio 
loudness does not depend strongly on M. Lacy et al. (2001) 



found no evidences for critical M or L/M ratio that turns on 
powerful radio jets. Ho (2002) showed that the dependence 
of radio loudness on M disappears altogether when one con- 
siders AGN with a broad range of intrinsic luminosity. At 
the same time he found that radio loudness seems to be re- 
lated to the mass accretion rate. RL and RQ sources are well 
separated in terms of M by the new Eigenvector analysis of 
Boroson (2002). 



The difference between RQ and RL on Eigenvectors di- 
agram (Boroson 2002) is probably connected with the inclu- 
sion of [Om]A5007 line parameters. In RL quasars part of 
the [Om] AA4959,5007 is coming from extended regions con- 
nected with the radio jet (Wilman, Johnstone & Crawford 
2000). Previous work has shown that RL activity in low M 
sources may be physically possible. Our sample contains a 
RL NLSyl that is one of the highest L/M sources, and other 
sources that radiate al log L/M « 4.5. Our analysis stressed 
that there are no appreciable effects (within the limits set 
by our S/N) on the H/3bc profile attributable to radio loud- 
ness, and that Pop. B RL and RQ sources can have the 
same M and L/M values. Therefore, it is not unreasonable 
to conclude that a similar range of M and L/M is physically 
possible for both RQ and RL sources. However, this does 
not mean that the mass function for RL and RQ sources 
is necessarily the same. A robust inference from the boot- 
strap analysis reported in H5.1l is that the mass function and 
the conditional probability of having certain L/M values at 
fixed M are likely to be different for the two AGN classes. A 
correct estimation of the mass function and of the L/M prob- 
ability distribution demands a more thorough analysis that 
includes selection biases. In addition, the RL/RQ dichotomy 
is probably related to parameters other than M and L/M, 
such as BH spin and morphological segregation. We suggest 
that the intrinsic mass function and L/M distribution dif- 
ferences as well as jet-related effects like [Om] AA4959,5007 
enhancement in RL sources and sample selection criteria 
may drive the RQ/RL separation shown by Laor (2000) and 
Boroson (2002). 
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Figure 8. H/3 line profiles after Fcllopt subtraction. From top 
to bottom the plots correspond to: (a) logL/M 4.4 - 4.9; (b) 
logL/M 3.9 - 4.4; (c) logL/M 3.4 - 3.9; (d) logL/M 2.9 - 3.4. 
Note that log L/M ^4.53 corresponds to L/LEdd ~ 1- Solid lines 
correspond to 7.3 $J logM^J 8.3, and dashed lines to 8.3 < logM^ 
9.3. 



8.3 Transition in H/3bc Line Profile Shape 

Connecting empirical and physical parameters, we find that 
median profiles for H/3bc in Pop. A sources (logL/M > 3.9, 
and for 3.5 < logL/M < 3.9 and 7.5< logL/M <8.5) show 
a Lorentzian shape. In other words, if logL/M >3.9, the 
shape is Lorentzian for all values of M. Similarly, if logL /M< 

3.4 (bins dominated by Pop. B surces), the profile is red- 
ward asymmetric and can be decomposed into two Gaussian 
components (one redshifted and the other not). This sug- 
gests that the transition from Lorentzian to double Gaus- 
sian profiles is governed by a critical L/M value (it is in- 
teresting to note that the H/fec asymmetry is also one of 
the major correlate of the original El of Boroson & Green 
1992). Unfortunately, the division of the bins in the L/M 
strip 3.4< logL/M <3.9 did not make the situation clearer. 
We conclude that the Lorentzian - double Gaussian pro- 
file transition occurs between 3.4< logL/M <3.9 (0.08 < 
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Figure 9. H/3bc line profile fitting after continuum, Fcllopt and 
H/3nc subtraction for the intervals in log L/M and M defined in 
Table 2. The "cleaned" H/3bc (dotted line) is shown along with 
single fitting components (dashed lines) and with the resulting 
best fit (thin solid line). 



L/L B dd <0.25). A more precise value would require knowl- 
edge of source orientation effects and a better-constrained a 
value. This provides a physical basis for the phenomenolog- 
ical finding about profile change at FHWM(H/3 B c)~ 4000 
km s _1 and the resulting Population A-B hypothesis (Sulen- 
tic et al. 2002). 



8.3.1 c(l/4) H/3bc Dependence on Mass: Not Only 
Gravitational Redshift 

We remark that a Double Gaussian fit to H/fec in Pop. B 
sources is a formal result: two components are required to ac- 
count for the redward asymmetry. The two components have 
a physical justification if, for example, the broader (VBLR) 
one can be ascribed to gas that lies closest to the continuum 
source in an optically thin (to the HI ionizing continuum) 
region with large covering factor f c ~ 1 (as defined in Su- 
lentic & Marziani 1993 and Brotherton 1996). A new result 
from this investigation is that the amplitude of the redward 
asymmetry is mass dependent. 

The c(l/4) dependence on M (for AlogM « 1 we ob- 



Searching For the Physical Drivers of Eigenvector 1 .-Influence of Black Hole Mass and Eddington Ratio 



- 1 - 

i„,„i.i,„ ?!5 
_ I I | I i i I | i i I i I i _ 


- 1 - 

; 1 : 

\ ;f i " ; 

if pi 

- 1 1 | 1 1 ! 1 | 1 1 1 1 "l" 1 - 


-_ - 




i ■ 





-5000 5000 
Av r [km s -1 ] 



-5000 5000 
Av r [km s" 1 ] 



Figure 10. Example of H/3bc line profile fitting with different 
functional forms. The upper panes show the H/3bc (dotted lines) 
fit by a Lorentzian (left, best fit identical to Fig.[5J and to a double 
Gaussian. The lower panel show residuals for both cases. The 
Lorentzian fit is significantly better, yielding to a decrease in x 2 
by a factor « 1.5. This is the most doubtful case among the ones 
considered in this study; in all other cases a Lorentzian/Gaussian 
exchange from best fits yields a \ 2 worsening by a factor 2. 
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Figure 12. Position of blue outliers (crosses) relative to other 
sources in our sample (open circles). The larger cross indicates 
the position of PKS 0736+01 - the only known RL "blue outlier" . 
The dashed line indicates the Eddington limit. Arrows indicate 
displacements of blue outliers if we apply an orientation correction 
~0.4 to the derived masses. Blue outliers are among the highest 
L/M sources and remain in the upper part of the diagram even 
after an orientation correction. One should consider that many 
other sources would move in the same direction of the blue outliers 
if a proper orientation were applied to all sources. 



tain a factor ~6 increase in the c(l/4) redward displace- 
ment) means that we can try to ascribe the presence of the 
"red shelf /VBLR component and the redward asymmetry 
to gravitational and transverse redshift (e. g., Corbin 1997). 
If the VBLR redshift is gravitational and transverse taking 
the c(l/4) value as a conservative estimate, we obtain the 
following distances from the central continuum source for 
the line emitting gas: if 3.5 < log L/M < 3.9, r « 0.005 pc 
and 0.01 pc for logM = 8 and logM = 9 respectively. In 
the case of the largest VBLR redshifts, the shift in radial 
velocity is Av r «2300 km s _1 and we find r w 0.015 pc 
~ 170 gravitational radii for logAf = 9. If we model the 
VBLR gas as a shell (f c ssl) with optical depth to the Ly- 
man continuum r < 1, a CLOUDY (Ferland 2000) simulation 
shows luminosity of H/3 logL(H/3)vBC ~ 41.7, where the line 
luminosity is in ergs s _1 . This falls far short in explaining 
the VBLR luminosity for sources in bin 3.5< log L/M < 3.9, 
8.5 < logM < 9.5. The average logL(H/3 B c) is 43.08. The 
VBLR contributes 2/3 of the total, so the VBLR average 
luminosity is logL(H/3)vBC ~ 42.9. The difference between 
the expected and observed VBLR luminosity is largely a 
consequence of the small shell radius required to explain 
the large Av r in the c(l/4). In addition, the assumption of 
c(l/4) being dominated by gravitational + transverse red- 
shift implies that z grav « 3/2(FWHM/c) 2 , which is not ob- 
served. We conclude that, even if c(l/4) is mass dependent, 
the c(l/4) shift amplitude cannot be explained by gravita- 
tional + transverse redshift alone. Once the stringent (and 
perhaps unphysical) distances from the central continuum 
sources set by gravitational redshift are relaxed, we note 
that our VBLR model can explain H/3 VBLR luminosity 
with reasonable shell radii. While this issue requires further 
investigation, it is tempting to suggest that we may be ob- 
serving the long-sought infall, and that we are observing it 
easily because of the low optical depth expected also along 
the line of sight to H/3. This would straightforwardly explain 
the ubiquitous presence of a strong H/3bc redward asymme- 
tries at low L/M. 



8.4 On the Nature of the Blue Outliers 

We observe "blue outliers" (following the simple wind model 
of Zamanov et al. 2002) when the [Om] AA4959,5007 lines 
and the peak of H/3 arise in different regions - H/3 from a 
near face-on accretion disk and [OIII] from a wind with ve- 
locity of about 1000 km s~ . In terms of orientation these 
RQ sources would be analogous to some CD RL quasars or 
even to BL Lac (blazars). Other nearly-as-large L/M ori- 
ented sources can/must exist in our sample but they are not 
"blue outliers" for at least two reasons: (i) they possess a well 
developed Narrow Line Region (NLR) which is unlikely to 
show a significant blue-shifted component; (ii) they are not 
oriented face-on. There are probably also lower L/M sources 
oriented face-on, however their L/M ratio may not be suffi- 
cient to power a wind or the wind outflow velocity may be 
low (roughly < 300 km s" 1 ). This could explain why the blue 
outliers are exclussively sources with very high L/M. It is 
worth noting that the only RL blue outlier, PKS 0736+01, 
is a flat spectrum radio source, and one of the brightest 
High Polarization, Optically Violently Variable sources. Its 
optical continuum is dominated by a "blazar" component 
as in BL Lac (Malkan & Moore 1986). Continuum proper- 
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ties suggest relativistic beaming making this a true pole-on 
RL source consistent with the interpretation of the "blue 
outliers" proposed by Zamanov et al. (2002). 

In the most recent scenarios for joint AGN and galaxy 
evolution (e.g., Granato et al. 2001), the active nucleus and 
the galaxy evolve together, with BH accreting matter and 
the galaxy making stars and supplying fuel for the quasar. 
At some point, the wind from the accreting BH blows away 
the matter surrounding it and a quasar emerges. The cen- 
tral accretion source then appears as an unobscured quasar 
which lasts as long as there is fuel in the accretion disk 
(Fabian 1999). In this scenario our "blue outliers" could 
represent the stage when the quasar has just emerged (cf. 
Krongold et al. 2001 for a similar interpretation stemming 
from the analysis of NLSyl host galaxies and environment). 
They could be in process of building of their NLR. 



9 CONCLUSIONS 

In this paper we computed virial masses and Eddington ra- 
tios for a sample of ~300 AGN. We have shown that L/M 
seems to govern the overall shape of the H/3bc, and we found 
an interesting dependence on M of the H/3bc asymmetry. 
The transition between sources showing Lorentzian and dou- 
ble Gaussian H/3bc profiles is likely associated with a critical 
L/M value, as suggested by Sulentic et al. (2002). L/M and 
M values in the range 8 < M < 9.5 and 3 < log L/M< 4.2 
are found likely for both RL and RQ AGN while others (log 
L/M> 4.2) may very unlikely for RL sources. However, it is 
important to stress the existence of a radio-loud blue out- 
lier radiating at very high L/M (log L/M~ 4.8; other RL 
sources radiate at logL/M ~4.4). This shows that there is 
no physical impossibility for RL sources to be of large L/M 
and small M; they may be simply less likely to be that way. 

We confirm that the "blue outliers" are mainly NLSyl 
sources (Zamanov et al. 2002; see also Sulentic et al. 2000a). 
We increase the number of known "blue outliers" to 12. We 
show that all blue outliers are accreting at 2 times higher Ed- 
dington ratio than other sources with similar mass. The com- 
pactness of their NLR points toward very young ages (Za- 
manov et al. 2002) , suggesting that they may be "fledgling" 
AGN as discussed for NLSyl by Sulentic et al. (2000a) and 
Mathur (2000). An interesting issue is then whether differ- 
ent M and L/M distributions may point toward evolution- 
ary effects, i.e., whether some RQ sources may be the parent 
population for all quasars (for example, NLSyls, "extreme 
Pop. A" sources, may evolve into Pop. B RQ or even RL 
depending on accretion of angular momentum and/or host 
galaxy type/evolution). 

In the future we will need of an orientation indicator for 
each individual AGN. This is very important for a correct 
evaluation of M. The CivA1549 line profile shows promise as 
an orientation indicator at least for Pop. A sources (Richards 
et al. 2002). It may also provide us with important clues 
about aspects of Broad Line Region structure and its de- 
pendence on L/M, and ultimately, with a 3D observational 
space to uniquely map into a 3D "physical" parameter space 
defined by M, L/M and orientation. 
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